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I. Nomenclature 


b - Unit vector in direction of local mag- 

netic field. 

9 P ~ Proportional control gain. 

9d - Derivative control gain. 

9i ~ Integral control gain. 

9w ~ Wheel speed control feedback gain. 

G - Gravity gradient torque. 

h - Spacecraft internal angular momentum. 

J - Spacecraft inertia matrix, 

n - Momentum management axis. Unit 

vector along intersection of magnetic 
control subspace and wheel control 
subspace. 

n± - II 2 space perpendicular to h. 

q - Attitude quarternion. 

Ri - Saturation ratio of ith actuator. Ri = 

|7v | /r max . 

v - Linear control of dynamics-inverted 

system. 

w - Wheel speed vector [ 0 w y w z ]'. 

v ) y , w z - y and z momentum wheel speeds, 

respectively. 

x - Unit vector along spacecraft x-axis. 

tp - Euler axis attitude error derived from 

e ? . e p = 2c ? (0)[ c ? (l) e q { 2) e 9 (3) ]'. 

tq - Quarternion error. = q^q^j. 

£w - Wheel speed error. e w = w — w re j. 

- Spacecraft angular velocity error. e w = 
w - u rc j. 

t - Spacecraft control torque. 

r i “ Attitude control torque projected onto 

h. 

r 2 * Attitude control torque projected onto 

n ± * 

Ti - Actuation command for ith actuator. 

Tmaxi “ Saturation value of ith actuator. 

T mi “ Attitude control torque projected onto 

effectuated by magnetic torquers. 
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Tm 2 ~ Momentum management torque paral- 

lel to h effectuated by magnetic dipole 
torquers. 

Tyj i - Attitude control torque projected onto 

nj. effected by momentum wheel 
actuators. 

r W 2 - Momentum management control torque 

parallel to n effectuated by momentum 
wheels. 

w - Spacecraft angular velocity. 

® - Quarternion multiplication operator. 


II. Introduction 

The Commercial Experimental Transporter 
(COMET) satellite provides a microgravity en- 
vironment for experimental payloads. For the 
first thirty days of the mission the satellite is op- 
erated in a solar-inertial fixed attitude to max- 
imize the power profile. The spacecraft is then 
pointed along the velocity vector and the exper- 
iment carrier is released for deboosting and re- 
covery. The spacecraft is then nadir pointed and 
additional microgravity experiments are con- 
ducted. 

The COMET Attitude Determination and 
Control System (ADACS) supports these mis- 
sion requirements by combining inverse dynam- 
ics nonlinear feedback with a novel momentum 
management technique. 

Inverse dynamics was developed in the mid- 
dle 1980’s for multi-mission, multi-configuration 
aircraft [1, 2] and subsequently applied to space- 
craft [3]. 

The technique uses nonlinear feedback to glob- 
ally transform nonlinear system dynamics into 
a linear system described by a linear second- 
order differential equation. Since the transfor- 
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mation is global and exact only a single set of 
control gains is required for both small-angle at- 
titude maintenance and large- angle attitude ma- 
neuvers. The technique inherently accomodates 
multiple spacecraft mass configurations without 
the need for multiple sets of control gains. The 
inverse dynamics technique compensates for and 
exploits spacecraft precession to achieve three- 
axis attitude control using two reaction wheels. 

The COMET momentum management sys- 
tem controls reaction wheel speed while pre- 
serving the satellite microgravity environment. 
Unlike the inner-loop/outer-loop structure com- 
mon to most satellites, the COMET momentum 
management system calculates the momentum 
management torques as an integral part of the 
attitude control torque. 


feedback using a Kalman observer to estimate 
the unmeasured velocity state. Including the 
position and velocity estimates from the Kalman 
observer, the control law described by equations 
(2) and (3) is modified 

r = J(q)v *f u> x [J(q)u + h] 4- G(q) (5) 

v = — g p iq — 9d^w — 9i I (6) 

Jo 

where the * indicates the observer-estimated 
value. 

The torque vector derived from the control 
law, equations (5) and (6), is applied to the 
COMET microgravity attitude actuators. 

IV. Torque Distribution & Momentum 
Management 


III. Inverse Dynamics Control 


The spacecraft dynamics are written 

J (q)u + u>* [J (q)u) + h]+ (?(?) — T - ( 1 ) 

s v 

precession grav grad 


The plant dynamics described by (1) are in- 
verted and globally linearized 


T 


V 


J{q)v + w* [J (q)uj + h] + G{q ) 


~9p € q 




9i / ff(A)dA. 

Jo 


( 2 ) 

(3) 


Substituting equations (2) and (3) into (1) 
yields the linear differential equation 


Cj = ~9p(q ~ 9dtw 



(4) 


The feedback gains, g p , gd, and gi are chosen 
using LQR design methods. 

Implementation of the inverse dynamics con- 
trol law described by equations (2) and (3), re- 
quires full-state feedback of spacecraft attitude 
and angular velocity. Many spacecraft lack the 
inertial measurement unit (IMU) necessary to 
measure angular rates. COMET has an IMU 
but its use is necessarily limited to conserve 
power. 

It was shown [4, 5] that inverse dynamics can 
be robustly implemented with only partial-state 


A. The Actuator Suite 

The COMET microgravity actuator suite con- 
sists of three 100 a-m ? electromagnetic torque 
coils along the spacecraft x, y, and z body 
axes and two 2 newton-meter-second momen- 
tum wheels along the spacecraft y and z body 
axes. There is no momentum wheel parallel 
to the spacecraft z-axis. (COMET is equipped 
with six 44 newton (10 lb) cold gas thrusters. 
They impart accelerations which exceed micro- 
gravity limits so their use is resticted to initial 
acquisition and recovery system reaction con- 
trol.) 

At the nominal 250 nautical mile altitude the 
magnetic torques deliver approximately 1 X 10 -3 
newton-meter torque in the plane perpendicu- 
lar to the local magnetic field. The momentum 
wheels deliver approximately 8.5 x 10 -3 newton- 
meter torque. 

B. Momentum Management by Torque Decom- 
position 

Three wheel and four wheel momentum man- 
agement schemes are usually implemented as a 
feedback outer-loop external to the attitude con- 
trol system inner-loop controller. Gas jets pro- 
vide instantaneous three-axis desatuation con- 
trol authority. For non-equatorial orbits mag- 
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netic dipoles provide instantaneous two-axis de- 
saturation and over-orbit three-axis momentum 
management. 

Inner-loop/outer- loop control structure sepa- 
ration works with with spacecraft having three 
or four wheels because the wheel axes span 
7Z 3 . The attitude control torque is exactly de- 
composed onto the wheel axes and effectuated 
entirely with the momentum wheel actuators. 
The instantaneous momentum error is projected 
onto the magnetic control subspace 7 l 2 normal 
to the local external field and is removed. Over 
orbit the magnet control subspace spans 7Z 3 
and wheel momentum is maintained. The con- 
trol problem separates whereby high-bandwidth 
attitude control is effected by the momentum 
wheels and low-bandwidth momentum control 
is effect by the magnetic torque coils. 

Having only two wheels, the COMET mi- 
crogravity control system cannot be partitioned 
into this inner-loop/outer-loop structure. The 
COMET attitude control torque must be ap- 
plied using a combination of the wheels and 
torque coils while momentum is simultaneously 
maintained over orbit. The momentum control 
torques are calculated continuously and as an 
integral part of the attitude control torques. 

The wheel torque subspace spans 'll 2 . When 
the spacecraft z-axis is not colinear with the lo- 
cal magnetic field the magnetic torque subspace 
spans a different 7 l 2 and the actuator torque 
space is overdetermined. The intersection of the 
two subspaces represents a direction in 7l 3 where 
the attitude control torque can be applied either 
using the magnetic torquers or the momentum 
wheels. This axis is the “momentum manage- 
ment axis” where some portion of the momen- 
tum error can be removed without disturbing 
spacecraft attitude. See figure 1. 

The COMET momentum management and 
attitude control torque distribution technique is 
captured mathematically by solving a system of 
linear constraints. The decomposition steps are 
summarized: 

£ x b 

n — — , 

\x x b\ 

r\ = nr. 

T2 = (/ — hh')r 


( 7 ) 

( 8 ) 
( 9 ) 


b 



n 


Figure 1: The “torque space” and the momen- 
tum management axis. 
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( 10 ) 


(ii) 


Equation (10) represents the constraints for solv- 
ing for the torque in subspace n ± . These constraints 
are: 

1. Zero magnetic torque in direction of b. 

2. Zero wheel torque in direction of x. 

3. The sum of the wheel and magnetic torque in 
their respective subspaces equals the compo- 
nent of the attitude control torque in space 
n x . 

4. Zero magnetic torque in direction of n. 

5. Zero wheel torque in direction of h. 

The 7x6 matrix in ( 10) is a rank 6 matrix and 
the equation has a unique solution. (In COMET 
these equations are solved using the numerically 
robust QR decomposition.) Equation (11) is un- 
derdetermined. Any combination of magnetic 
and wheel torque in the direction h which equals 
ri will correctly maintain the spacecraft attitude. 

The momentum management torque is pro- 
portional to the wheel speed error projected 
onto h. A wheel speed corrective torque and 
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a compensating magnetic torque are applied in 
the n and -n direction, respectively. 

t w 2 = —g w nh't w ( 12 ) 

r m2 = T \ ~ T w 2 - ( 13 ) 

The final actuator torques are the sum of the 
components in the h and ni spaces 

T w — T w 1 + T w 2 (I'l) 

T m = rVni T m 2- ( 1 ^) 


C. Saturation Control 

The computed wheel torques and magnetic 
dipole moments may drive the actuators into 
saturation. Saturation does not usually occur 
uniformly across all the actuators. The effect of 
saturation is to misdirect the corrective rotation 
of the satellite. 

Saturation control in COMET is based upon 
the heuristic that it is preferrable to maintain 
the direction of the control torque at reduced 
magnitude than to apply maximum control au- 
thority in the wrong direction. 

The five actuator commands (two wheel 
torques, three magnetic dipole values) are calcu- 
lated and five saturation ratios il, = |T,|/T mar i 
are calculated. If any of these five ratios exceeds 
unity than all the actuator commands are scaled 
by the largest of these five ratios. 

The effect of occasional actuator saturation 
upon linear systems is generally benign and not 
destabilizing. The effect of actuator saturation 
upon an inverse dynamics control system is com- 
plicated in that saturation leads to incomplete 
inversion. Simulation has demonstrated good 
system performance even with occasional satu- 
ration. The effect of saturation on incomplete 
dynamics inversion and control system perfor- 
mance is an area which requires future research. 


V. Simulation 

Simulations demonstrate the performance 
of the COMET microgravity control system. 


Gravity gradient and aerodynamic torques are 
the dominant distubances and are simulated. 
The simulation orbit is circular at 463 km (250 
NM) altitude. The F10.7 solar flux level is 116. 

Spacecraft state estimates generated by atti- 
tude determination system are used for calcula- 
tion the attitude and rate errors. Except where 
noted, the IMU is off and feedback linearization 
is achieved using angular velocity estimates. 

Simulation series (1) demonstrates COMET 
solar-pointing performance with a momentum 
bias. The spacecraft attitude recovers from the 
initial 15° error. Momentum wheel rates con- 
verge to their nominal 1500 RPM bias levels. 

Lack of sensor input combined with unmod- 
elled aerodynamic torques cause the attitude er- 
ror to increases in the umbra. This pointing er- 
ror is due to attitude estimation error developed 
in the Kalman Filter extrapolation step. Simu- 
lation series (2) shows represents the same run 
as series (1) except the IMU is on and satellite 
body rates are directly measured. Umbra per- 
formance is greatly improved. 

Simulation series (3) demonstrates COMET 
solar-pointing performance without momentum 
bias. The spacecraft attitude recovers from the 
initial error. Momentum wheel mean rates con- 
verge to zero. The attitude control system in 
series (3) is identical to that of series (1) demon- 
strating control system flexibility in accomodat- 
ing different momentum configurations. This 
flexibility is a consequence of the exact global lin- 
earization provided by the inverse dynamics con- 
trol law. With equal flexibility the control sys- 
tem accomodates spacecraft with varying mass 
properties. COMET is a multi-configuration 
spacecraft which uses the same control law for 
all its mass configurations. 

Simulation series (4) demonstrates COMET 
Earth-pointing performance 1 . 

Simulation series (5) demonstrates 
the COMET large-slew microgravity transition 

'In Earth-pointing mode the COMET spacecraft coor- 
dinate system differs from the conventional local verti- 
cal system. In this mode the roll axis (x-axis) is nadir- 
pointing, the yaw axis (x-axis) is aligned anti-parallel 
to the orbit normal, and the pitch axis (y-axis) com- 
pletes the triad, pointing approximately in the velocity 
direction. 
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from solar-pointing to Earth-pointing. Note the 
rate quenching of the COMET pitch momentum 
wheel (yaw wheel by conventional terminology). 
Spacecraft acceleration is maintained within the 
10 _5 jf microgravity acceleration limit during the 
slew. 

Simulation series (6) demonstrates COMET 
slewing 90° in pitch and yaw and holding to an 
arbitrary attitude. Momentum wheels maintain 
their 1500 RPM bias during the slew. This mode 
is used for reentry attitude control prior to re- 
covery capsule deboosting. 

Simulation series (7) demonstrates the synthe- 
sized roll control performed without a roll wheel. 
This synthesized roll is a consequence of the in- 
clusion of the precession term in the inverse dy- 
namics. Series (7a) shows the large roll when 
the precession term is not included. Simulation 
series (7b) shows the roll error reduction caused 
by inclusion of the precession term. 

Simulations (l)-(7) demonstrate the perfor- 
mance of the COMET attitude control sys- 
tem. Slewing of the the reaction wheel during 
transition and maintenance of their speeds in 
the steady-state demonstrates the performance 
of the novel COMET momentum management 
technique. Operational flexibility to the va- 
riety mass properties, momentum biases and 
command attitudes is a feature inherent to the 
method of inverse dynamics control. 

From simulation experiments, inverse dynam- 
ics control is empirically known to be robust to 
large variations in the “model” plant. Simu- 
lation series (8) demonstrates that COMET is 
robust to large uncertainty to mass properties. 
In this simulation the three mass principle axes 
of the true satellite are skewed 15° from their 
“model” orientation. The system remains sta- 
ble and within specification. The mass proper- 
ties mismatch does causes a 50% increase RMS 
attitude error (see simulation series 1 for com- 
parison.) 


VI. Conclusion 

The COMET attitude determination and con- 
trol system, using inverse dynamics and a 


novel torque distribution/momentum manage- 
ment technique, has shown great flexibility, per- 
formance, and robustness. 

Three-axis control with two wheels is an in- 
herent consequence of inverse dynamics con- 
trol which allows for reduction in spacecraft 
weight and cost, or alternatively, provides a sim- 
ple means of failure-redundancy for three-wheel 
spacecraft. 

The control system, without modification, has 
continued to perform well in spite of large 
changes in spacecraft mass properties and mis- 
sion orbit altitude that have occurred during 
development 2 . This flexibility has obviated im- 
position of early stringent ADACS design con- 
straints and has greatly reduced commonly in- 
curred ADACS modification costs and delay as- 
sociated with program maturation. 
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increases in spacecraft mass has lead to reduction of 
orbit altitude from the original 556 km (300 nm) al- 
titude to the 463 km (250 nm) altitude. The control 
system was designed for the original altitude. 
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